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SUMMARY 

The  Lennard- Jones  6-12  atom  interaction  potential  is  applied  to  the  cal¬ 
culation  of  adsorption  energy  of  an  atom  onto  as  many  as  132  sites  on  a  unit 
cell  area  of  each  of  the  eight  highest  surface  density  planes  of  a  body- 
centered-cubic  substrate.  From  the  calculations,  topographical  maps  of  adsorp¬ 
tion  energy  on  a  unit  cell  area  of  each  plane  for  various  adsorbate-adsorbent 
combinations  are  plotted.  Normalized  values  of  maximum  adsorption  energies  and 
minimum  surface  diffusion  activation  energies  for  all  cases  are  tabulated.  The 
results  are  compared  with  available  experimental  data  concerning  adsorption  of 
alkali  metals,  alkali  earth  metals,  and  inert  gases  on  transition  metal  sub¬ 
strates  and  are  in  good  agreement.  Application  of  the  results  to  transition- 
metal  -  transition-meoal  combination  predicts  that  the  110  surface  is  the 
lowest  energy  configuration  for  a  body-centered-cubic  crystal. 


INTRODUCTION 


Adsorption  on  an  atomic  scale  gained  recognition  with  the  advent  of  the 
field-ion  microscope  and  field-emission  microscope,  and  it  gained  in  importance 
proportionally  with  the  fields  of  ultrahigh  vacuum,  gas -filled  thermionic  con¬ 
verters  and  the  contact-ionization  engine.  Through  the  use  of  field-ion  and 
field-emission  microscopy,  it  has  been  observed  that  adsorption  properties  are 
strongly  dependent  on  the  atomic  arrangement  of  the  substrate  (in  some  cases 
more  so  than  on  the  materials  involved) .  Forces  governing  adsorption  of  many 
species  on  metals  or  semiconductors  have  been  found  ta  be  of  the  form  of  dis¬ 
persion  forces  such  as  van  der  Waal  forces  /(ref,  l).  /  By  considering  dispersion 
force  interactions  of  an  adsorbing  atom  with  the  nearest  neighbors  of  the  sub¬ 
strate,  one  can  verify  that  the  substrate  structure  grossly  affects  the  adsorp¬ 
tion  propertiesT(refs.  2  and  3).  / It  is  the  intent  of  this  report  to  calculate 
the  interacbrCn  potential  encountered  by  a  single  atom  at  any  point  on  various 
ideal  single  crystal  planes  of  the  body-centered-cubic  structure.  Furthermore, 
it  is  intended  to  deduce  from  such  potential  surfaces  some  general  properties 
of  physical  adsorption  and  specific  properties,  such  as  adsorption  energy  and 
diffusion  activation  energy,  and  to  compare  these  results  with  experimental  ob¬ 
servations^  _ _ 


frhe  Lennard-Jones  6-12  potential  is  used  as  a  model.  Its  applicability 
to  adsorption  of  a  nonpolar^  atom  on  a  metal  is  discussed  in  the  following  sec¬ 
tion.  The  general  mathematical  techniques  employed  are  outlined. in  the  section 
PROCEDURE,  while  details  concerning  the  performance  of  the  specific  lattice  sum¬ 
mations  for  the  body-centered-cubic  structure  are  contained  in  appendix  A.  The 
last  section  compares  the  results  with  the  scant  experimental  data  and  discusses 
the  application  of  the  results  to  known  adsorption  characteristics  of  practical 
surfaces 
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MODEL 


The  Lennard-Jones  potential  is  (ref.  2) 


E  = 


where  E  is  the  atom-atom  interaction  energy,  e  is  the  depth  of  the  energy 
well,  a  is  that  finite  value  of  r  for  which  E  is  zero,  and  r  is  the  in- 
ternuclear  distance.  The  attractive  portion  of  equation  (l)  is  proportional  to 
the  inverse  sixth  power  and  represents  a  long-range  van  der  Waals  interaction, 
while  the  repulsive  contribution  is  approximated  by  an  inverse  twelfth  power  de¬ 
pendence. 


Application  of  this  atom-atom  interaction  potential  to  adsorption  on  a 
solid  assumes  additivity  of  the  interaction  of  the  adsorbate  with  each  atom  of 
the  adsorbent  or 

9  ’  4€Z  tfe)6  +  (e)1  (2) 

i 

where  cp  is  the  atom-metal  interaction  energy  and  r^  is  the  distance  from 
the  adsorbed  atom  the  i^*1  atom  of  the  substrate.  Thus,  effects  arising  specif¬ 
ically  from  the  structure  of  the  substrate  are  contained  in  the  lattice  summa¬ 
tions,  rT®  and  ^  r7-^  whereas  effects  arising  from  the  atomic  properties 
i  i 

of  the  constituent  atoms  are  contained  in  the  parameters  a  and  e.  Appen¬ 
dix  B  contains  a  discussion  of  the  physical  meanings  of  both  a  and  e. 

Equation  (2)  neglects  any  contribution  to  the  interaction  by  free  electrons 
in  the  solid;  however,  a  study  by  Pierotti  and  Halsey  (ref.  4)  on  the  inter¬ 
action  of  krypton  atoms  with  metals  has  shown  that  theories  which  include  the 
free  electron  contribution  of  the  metal  do  not  give  significantly  improved  re¬ 
sults  as  opposed  to  theories  in  which  the  free  electron  contributions  are  ne- 


-*-All  atoms  have  instantaneous  dipole  moments  that  are  the  primary  source 
for  the  long  range  interactions  between  atoms.  An  atom  is  considered  nonpolar 
if  when  isolated  its  net  dipole  moment  is  zero  over  a  period  of  time. 
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glee ted.  Neustadter,  Luke,  and  Sheahan  (ref.  5)  have  further  shown  that  very- 
low  coverage  adsorption  of  various  alkali  raetal  atoms  on  tungsten  can  be  satis¬ 
factorily  explained  by  using  equation  (2).  Of  course,  the  validity  of  the  pro¬ 
posed  model  can  best  be  tested  by  comparing  the  calculated  results  with  experi¬ 
mental  data. 


PROCEDURE 


In  order  to  preserve  the  general  applicability  of  this  work  to  adsorption 
on  any  substrate,  equation  (2)  is  rewritten  as 


-i-  = 

4e 


(3) 


where  d-  is  r*/a.0,  and  aQ  is  the  lattice  constant  of  the  substrate.  Equa¬ 
tion  (3)1is  treated  as  follows:  For  each  surface  orientation  considered,  a  co¬ 
ordinate  system  is  chosen  such  that  the  x,y  plane  passes  through  the  centers 
of  all  the  surface  atoms  and  z  is  positive  along  the  outward  normal.  Clearly 
<p/4e  is  a  function  of  x,  y,  and  z. 


As  shown  in  sketch  (c)  of  appendix  A  (p.  9)  each  unit  cell  is  broken  up 
into  a  grid  of  as  many  as  132  prospective  adsorption  sites  where  the  number  of 
sites  is  chosen  so  that  each  unit  cell  will  have  approximately  the  same  density 
of  prospective  sites.  At  any  one  of  these  sites  (x^,y^)  the  sums  dp  and 


(b) 


i 

]£^d^2  are  evaluated  as  functions 
i 

of  z.  Using  these  sums  in  equa¬ 
tions  (3)  and  plotting  cp/4e  as  a 
function  of  z  yield,  for  example, 
sketch  (a)  (which  will  depend  on  sur¬ 
face  orientation,  site,  and  a/a0 
for  its  specific  form). 

The  assertion  is  made  that,  if 
adsorption  were  to  occur  over  site 
xi,yi,  the  nucleus  of  the  particle 

would  be  located  at  zm±n,  and  the 
bond  energy  would  be  Tmin/46  (where 
the  quantities  are  as  defined  in 
sketch  (a)).  The  derivative  of  cp / 4€ 
is  taken  (numerically)  with  respect 
to  z ,  and  cpinin/ ^  and  ^min  are 
calculated  for  a  specific  c/a0  (see 
appendix  A).  This  entire  process  is 
now  repeated  for  each  predetermined 
location  Xj,yj  on  a  unit  cell  of 
the  adsorbent.  Then  9min/4e  is 
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plotted  as  a  function  of  surface  coordinates  x  and  y.  An  example  of  such  a 
plot  is  shown  in  sketch  (b)  for  the  110  plane  and  o/a0  =1.00  (corresponding 
to  barium  on  tungsten).  Here  one  unit  cell  of  surface  area  is  represented  with 
the  locations  indicated  for  the  substrate  atoms  nearest  the  surface.  Along 
each  of  five  equally  spaced  lines  of  constant  y,cpmln/4e  is  plotted  against  x 
with  the  ordinate  cpmin/ 4e  in  the  y-direction. 

The  minimum  value  of  cpmin/ 4e  is  read  graphically  and  is  listed  in 
table  I.  These  values  are  the  normalized  heat  of  adsorption  $/4e  for  the  pre¬ 
ferred  adsorption  site,  for  a  specific  orientation,  and  for  a  specific  value 
of  cr/a0.  The  quantity  Eq/ 4e  is  also  determined  graphically  by  identifying 
it  with  the  least  amount  of  energy  that  must  be  supplied  to  an  adsorbed  atom  to 
enable  it  to  move  from  one  preferred  adsorption  site  to  another,  where  the 
assumed  path  was  not  necessarily  a  straight  line.  In  the  110,  210,  and  310 
cases,  the  topographical  potential  plots  were  obtained  by  manually  cross - 
plotting  the  data  for  cpmin/ against  position.  In  the  remaining  cases,  the 
topographical  potential  plots  were  generated  directly  by  a  Fourier  interpolation 
7090  program  from  the  original  digital  output.  This  procedure  was  made  possi¬ 
ble  by  using  the  fact  that  the  potential  is  periodic  along  any  line  parallel  to 
a  border  of  the  unit  surface  cell.  This  entire  procedure  is  repeated  for  a 
range  of  values  of  a/aQ  and  the  remaining  substrate  orientations.  For  the 
specific  summation  technique  applied  to  the  body-centered-cubic  structure  see 
appendix  A. 


DISCUSSION  OF  RESULTS 

The  minimum  values  of  cp/4e  and  Eq/4e  are  listed  in  table  I.  The  orig¬ 
inal  plots  of  the  potential  variation,  such  as  illustrated  in  sketch  (b),  were 
replotted  as  topographical  maps  showing  lines  of  equal  adsorption  energy.  The 
eight  highest  surface  density  planes  were  considered  (110,  100,  211,  310,  111, 
321,  411,  and  210  in  order  of  decreasing  density)  and  are  presented  in  fig¬ 
ures  1  to  8,  respectively.  Included  in  the  topographical  maps  are  the  locations 
of  the  substrate  atoms  in  each  unit  cell.  The  most  probable  migration  direction 
is  assumed  to  be  path  of  least  resistance  (minimum  Eq/4e)  connecting  two  equiv¬ 
alent  adsorption  sites. 

The  values  of  a/aQ  used  in  the  calculations  were  0.80,  0.894,  1.00, 

1.123,  and  1.20.  Table  II  lists  the  values  of  cr/a0  corresponding  to  various 
adsorbent-adsorbate  combinations.  The  most  meaningful  interpretation  of  these 
data  results  from  considering  the  comparative  effects  of  substrate  surface 
structure  (crystal  orientation)  on  adsorption  and  diffusion  for  a  given  combi¬ 
nation  of  adsorbent  and  adsorbate  (i.e.,  for  a  given  a/ aQ) . 

The  Lennard-Jones  potential  model  corresponds  most  closely  to  inert  gas 
atoms  adsorbing  on  a  semiconductor  or  an  insulator  when  the  coverage  is  low 
enough  to  eliminate  any  interaction  of  adsorbate  atoms  with  one  another.  This 
limit  of  coverage  is  shown  experimentally  in  some  cases  to  be  as  high  as  0.7 
(e.g.,  thorium,  ref.  6),  thus  allowing  the  calculations  to  be  useful  in 
practical  systems  at  coverages  other  than  zero.  The  Lennard-Jones  potential 
model  has  been  used  in  calculating  adsorption  energies  of  various  metals  on 
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tungsten  (ref.  7)  and,  in  particular,  adsorption  energies  of  alkali  metals  on 
tungsten  (ref.  5) .  Agreement  of  the  results  with  experimental  data  implies 
that  this  method  may  also  be  applied  to  other  metal -metal  systems.  Substrates 
that  are  compared  are  assumed  to  be  atomically  smooth  ideal  crystal  planes  with 
no  dislocations,  steps,  impurities,  missing  atoms,  or  surface  relaxations.  The 
experimental  situation  most  closely  approaching  these  assumptions  exists  in  the 
field-emission  microscope  in  which  the  tip  is  an  annealed  single  crystal  with  a 
minimum  of  imperfections  and  on  which  there  are  small  areas  of  atomically  smooth 
planes . 

In  emission-microscope  experiments  by  Utsugi  and  Gomer  (ref.  8)  and 
Drechsler  (ref.  9)  who  studied  the  adsorption  of  barium  on  tungsten,  the  acti¬ 
vation  energies  for  diffusion  Ed  on  the  110  and  321  planes  were  measured  and 
were  in  good  agreement  with  the  calculated  results  as  follows: 

Measured: 


(11Q)Ed  0.41 

(321)Ed  0.83 


=  0.49  (Ref.  8) 


(llO)Ed  0.20 
1321 )Ed  “  0.65 


0.31  (Ref.  9) 


Calculated: 

(llO)Ed/4€ 

(32l)Ed/4e 

Experimental  data  such  as  these  may  be  used  to  evaluate  e  empirically,  which 
in  turn  may  be  used  to  obtain  quantitative  values  for  Ed  and  $  for  each 
plane  of  a  given  system  (i.e.,  given  a/aQ) .  An  alternative  theoretical  expres¬ 
sion  for  e  is  presented  in  appendix  B. 

Calculated  adsorption  energies  from  table  I  (with  a/a^  =  1.00  for  xenon 
on  tungsten  and  molybdenum)  agree  favorably  with  emission  microscope  experi¬ 
ments  by  Ehrlich  and  Hudda  (ref.  10) ,  who  studied  the  adsorption  of  xenon  on 
tungsten  and  molybdenum.  They  found  stability  of  the  adsorbate  highest  on  the 
210,  310,  and  611  planes,  lower  on  the  211  plane,  and  lowest  on  the  100  and  110 
planes  at  temperatures  where  diffusion  is  not  important.  This  remarkable  agree¬ 
ment  is  less  surprising  when  one  considers  that  the  model  is  best  suited  to 
inert  gas  interactions. 

Adsorption  of  cesium  on  transition  metals  is  of  considerable  interest  in 
those  applications  related  to  electron  and  ion  emission.  Unfortunately,  the 
adsorption  energy  of  cesium  on  various  planes  of  metals  is  not  as  well  docu¬ 
mented  as  that  of  inert  gases  on  metals.  It  has  been  repeatedly  observed  (e.g., 
see  ref.  ll)  that  the  110  plane  of  both  molybdenum  and  tungsten  has  a  higher 
current  density  than  any  other  plane  at  a  given  substrate  temperature  and  ce¬ 
sium  arrival  rate  (in  the  region  of  low  to  intermediate  coverages) .  Since  this 
plane  has  the  highest  vacuum  work  function,  these  observations  have  been  inter¬ 
preted  as  implying  a  maximum  stability  of  cesium  on  the  110  surface,  whereas 


=  =  0.39 

0.345 
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the  calculations  predict  that  cesium  should  be  least  stable  on  the  110  plane  of 
tungsten  (a/a  =  1.123).  However,  one  may  attribute  increased  emission  density 
from  a  partly  covered  plane  to  any  one  or  combination  of  the  following  three 
factors : 

(1)  High  stability  of  the  adsorbate  on  the  plane  in  question,  which  de¬ 

pends  exponentially  on  $  and  E^ 

(2)  High  number  of  sites  available  for  adsorption 

(3)  Large  effective  dipole  moment  of  an  adsorbed  particle 

The  110  plane  has  the  highest  density  of  sites  available  of  all  body- 
centered-cubic  planes.  Good  and  Muller  (ref.  6)  have  shown  that  for  a  given 
surface  density  of  adatoms  the  110  showed  by  far  the  greatest  decrease  in  work 
function  of  all  body-centered-cubic  planes.  The  relatively  large  emission  den¬ 
sity  from  the  110  plane  of  tungsten  partly  covered  with  cesium  is  not  suffi¬ 
ciently  definitive  of  adsorption  energies  to  provide  an  adequate  test  of  the 
Lennard- Jones  potential  model. 

In  general,  the  calculations  substantiate  the  prevalent  model  of  surface 
adsorption  in  which  discrete  adsorption  sites  exist  and  in  which  diffusion  is 
defined  as  the  motion  of  atoms  between  well-defined  sites.  Also  the  calcula¬ 
tions  verify  the  expectation  that  these  sites  exist  at  wells  in  the  surface 
where  the  outer  electrons  of  adsorbate  participate  in  binding  with  a  maximum 
number  of  substrate  atoms. 

Additional  verification  of  the  validity  of  the  calculations  is  realized  in 
examining  self-adsorption  and  self-diffusion.  The  lowest  energy  surface  config¬ 
uration  of  a  body-centered-cubic  crystal  of  tungsten  (refs.  12  and  13)  was  shown 
experimentally  to  be  the  110  surface,  where  the  experimental  situation  closely 
approaches  that  assumed  in  the  calculations.  Recognizing  the  relatively  low 
values  of  $/4e  and  Ej/4e  at  a/aQ  =  0.80  (a/aQ  =  0.786  for  self-adsorption) 
for  the  110,  100,  and  211  planes  (in  that  order)  one  can  see  that  the  calcula¬ 
tions  predict  that,  for  a  body-centered-cubic  crystal  annealed  in  vacuum,  the 
final  surface  would  consist  of  110  planes,  while,  in  the  transition,  the  100 
and  211  planes  have  a  high  probability  of  existing.  It  must  be  noted  here  that 
even  though  the  physical  mechanism  of  metallic  binding  may  not  be  described  ex¬ 
actly  by  the  Lennard-Jones  potential,  the  structure-dependent  adsorption  prop¬ 
erties  of  such  a  system  can  be  readily  predicted  by  using  the  Lennard-Jones  po¬ 
tential  summation. 


SUMMARY  OF  RESULTS 

\  The  following  four  factors  serve  to  support  the  applicability  of  the  results 
to  practical  adsorption  phenomenon: 

1.  The  general  picture  of  adsorption,  as  observed  experimentally,  is  ver¬ 
ified  because  atomically  rough  surfaces  tend  to  have  higher  heats  of  adsorption 
than  closer  packed  surfaces.! 


rr  More  specifically,  the  order  of  heats  of  adsorption  for  xenon  on  var¬ 
ious  planes  of  tungsten  is  the  same  as  that  observed  in  field-emission- 
microscope  experiments. 

3.  The  equilibrium  surface  of  body-centered  cubic  crystals  is  the  110 
plane  determined  both  by  experiment  and  calculation. 


4.  The  relative  diffusion  activation  energy  for  barium  on  tungsten  is  in 
good  agreement  with  measured  values.  The  possible  extension  of  these  calcula¬ 
tions  to  the  adsorption  properties  of  other  structures  (e.g.,  face-centered- 
cubic,  diamond,  etc.)  and  of  single -crystal  surfaces  with  various  imperfections 
(e.g.,  missing  atoms  and  steps)  can  readily  be  seen^J 


Lewis  Research  Center, 

National  Aeronautics  and  Space  Administration, 
Cleveland,  Ohio,  August  6,  1965. 
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APPENDIX  A 


DETAILS  OF  THE  SUMMATION  FOR  BODY -CENTERED -CUBIC  SUBSTRATES 


The  location  of  all  atoms  in  a  body-centered-cubic  crystal  with  respect  to 
the  coordinate  system  described  previously  has  been  recorded  by  Bacigalupi 
(ref.  14)  for  twenty  orientations.  Thi£  technique  can  be  adapted  to  the  present 
problem  to  define  the  four  quantities  S,  X,  V,  and  d,  where  S  locates  the 
nearest  surface  atom  to  any  surface  atom,  X  locates  the  nearest  surface  neigh¬ 


bors  in  the  x-direction  to  any  surface 
atom,  V  connects  any  surface  atom  to 
the  closest  atom  in  the  next-to-surface 
plane,  and  d  is  the  interplanar  dis¬ 
tance  measured  along  the  z-axis.  The 
vectors  X  and  S  form  a  unit  cell 
of  the  surface  as  shown  in  sketch  (c). 

With  these  definitions  and  equa¬ 
tion  (Al)  all  the  atoms  in  the  crystal 
with  respect  to  the  surface  plane  can 
be  located 


Xj  =  xc  +  X±| V(X) J  +  n|X|  +  C  ' 

h  =  *o  +  hpMl  +”|sh)|  l  (ai) 

c  =  [yi-Yo-i1|t(x)|])||}^ 

where  X^Y^  is  the  location  of  any  atoms  referred  to  the  surface  plane  (x,y 
plane),  X0,Y0  is  the  starting  point  and  location  of  the  upper  left  corner  of 
the  unit  cell,  1^  is  the  plane  in  which  the  i^*1  atom  exists  (0  for  surface 
plane;  1  for  next-to-surface  plane,  etc.)  (x)  or  (Y)  indicates  the  component  of 
the  particular  vector,  and  n  and  m  are  integers. 


In  performing  the  summations,  it  was  necessary  to  establish  a  cutoff  for 
the  computer.  Using  an  integral  approximation  to  determine  the  contribution  of 
those  terms  omitted  from  the  summation  leads  to  the  result  that  accuracy  to 
within  1  percent  is  achieved  when  the  summation  is  extended  over  all  the  atoms 
within  seven  aQ  of  the  site  being  considered.  Therefore  equation  (A2)  is  in¬ 
troduced  to  establish  the  cutoff: 


XQ  +  7.0  +  |XJ  >  X±  >  XQ  -  7.0 


Yq  +  7.0  >  Y±  >  Y0  -  7.0  +  j  S(Y) j 


(A2) 


The  prospective  adsorption  sites  for  which  the  summations  were  performed 
are  located  at  the  intersections  of  five  to  ten  equidistant  lines  parallel 
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to  X  and  eleven  to  twenty  equidistant  lines  parallel  to  S,  as  shown  in 
sketch  (c).  The  distance  r^  from  an  adsorbing  atom  to  the  fth  atom  in  the 
crystal  can  now  be  evaluated  as  a  function  of  z.  At  each  site,  d^  (defined 
by  eq.  (A3)) 

ai  =  [fri  -  *oY  +  (y±  -  Y0)2  +  (lid  +  Z)2]l/2  (A3) 

was  evaluated  for  each  of  sixty  evenly  spaced  values  of  z  ranging  from  0.05 
aQ  to  2.1  aQ.  The  corresponding  values  of  dj_  were  used  in  the  summation  in 
equation  (3).  The  minimum  value  of  cp/4e  obtained  from  equation  (3)  for  a 
given  site  was  used  along  with  one  value  of  cp/4e  in  the  negative  z-direction 
and  three  values  in  the  positive  z-direction.  A  fourth-order  equation  was 
fitted  to  these  five  points,  and  a  minimum  was  determined.  The  minimum  was 
further  corrected  by  using  the  Newton-Raphson  method.  Sketch  (b)  is  an  ex¬ 
ample  of  these  minima  plotted  as  a  function  of  their  position  for  the  110 
plane  for  a/ a0  =  1.00. 


9 


APPEND  IX  B 


EVALUATION  OF  a/aQ  AND  e 

The  Lennard-Jones  potential  (ref.  2)  for  interactions  between  two  atoms  is 


E  =  4e 


[■  at  •  ef] 


(1) 


which  graphically  is  of  the  form 
shown  in  sketch  (d) ,  where  a  is 
that  finite  value  of  r  for  which 
E  =  0,  e  is  the  interaction  bond 
energy  of  the  two  atoms ,  and  rQ 
is  their  internuclear  equilibrium 
distance . 

To  evaluate  a  note  that  rQ 

is  that  value  of  r  for  which  E 

is  a  minimum.  Therefore ,  a  can 

be  related  to  rQ  by  setting 

(5E/dr)r  =  0.  This  gives 
x  o 

a  =  r02"l/6.  Since  r  can  be  ap¬ 
proximated  as  the  sum  of  the  hard- 
sphere  radii  of  the  two  interacting 


atoms,  values  can  now  be  calculated  for  o.  Table  II  lists  a/aQ  for  various 
combinations  of  adsorbates  and  adsorbents  with  body-centered-cubic  structures. 


Various  theories  have  been  developed  to  evaluate  explicitly  the  attractive 
long-range  interaction  energy  for  two  dissimilar  atoms.  The  best  results  seem 
to  be  obtained  (ref.  4)  from  the  theory  developed  by  Kirkwood  (ref.  15)  and 
Muller  (ref.  16) .  They  obtained 


E  (attractive) 


6mc2  aAaM 
r6  “A 
XA 


(Bl) 


where  a  is  the  electronic  polarizability,  X  is  the  diamagnetic  susceptibility, 
m  is  the  electron  rest  mass,  c  is  the  speed  of  light  in  vacuum,  the  subscript 
A  refers  to  the  adsorbed  atom,  and  the  subscript  M  refers  to  the  substrate 

metal  atoms . 

Comparison  of  equation  (Bl)  with  the  attractive  term  in  the  Lennard-Jones 
potential  allows  evaluation  of  e  from 


4ea^ 


6mc^ 


aAaM 


+ 

XA  % 


(B2) 
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2.79 
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2.719 

0.262 

3.009 

0.12 

3.35 

0.340 

2.772 

0.17 

3.04 

0.275 

3.417 

0.640 

3.071 

TABLE  II.  -  VALUES  OF  a/aQ  FOR  VARIOUS  ADSORBATES 
OH  BODY-CENTERED  CUBIC  SUBSTRATES a 


Adsorbate 

Substrate 

Vanadium 

Chromium 

Iron 

Tantalum 

Niobium 

Molybdenum 

Tungsten 

Neon 

0.864 

0.890 

0.891 

0.825 

0.828 

0.846 

0.844 

Argon 

0.955 

0.985 

0.987 

0.909 

0.911 

0.934 

0.932 

Krypton 

0.981 

1.013 

1.015 

0.933 

0.935 

0.960 

0.957 

Xenon 

1.040 

1.075 

1.077 

0.987 

0.990 

1.017 

1.013 

Lithium 

0.849 

0.874 

0.875 

0.812 

0.813 

0.832 

0.830 

Sodium 

0.951 

0.981 

0.982 

0.905 

0.908 

0.931 

0.927 

Potass ium 

1.084 

1.122 

1.124 

1.027 

1.029 

1.064 

1.061 

Rubidium 

1.227 

1.162 

1.164 

1.063 

1.066 

1.096 

1.092 

Cesium 

1.178 

1.220 

1.224 

1.114 

1.117 

1.150 

1.123 

Thorium 

0.921 

0.950 

0.951 

0.878 

0.880 

0.902 

0.894 

Barium 

1.045 

1.080 

1.083 

0.992 

0.994 

1.021 

1.018 

Strontium 

1.017 

1.050 

1.053 

0.966 

0.968 

0.994 

0.991 

Self¬ 

adsorption 

0.784 

0.785 

0.784 

0.786 

0.787 

0.787 

0.785 

aBased  on  atomic  radii  obtained  from  ref.  17,  p..  403. 
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Figure  2.  -  Topographical  map  of  normal¬ 
ized  adsorption  energy  on  100  plane. 
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(c)  o/aQ  =  1. 20. 


Figure  4.  -  Topographical  map  of  normalized  adsorption  energy  on  310  plane. 
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Figure  8.  -  Topographical  map  of  normalized  adsorption  energy  on  210  plane. 
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